Mott insulating and glassy phases of polaritons in ID arrays of coupled cavities 
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By means of analytical and numerical methods we analyze the phase diagram of polaritons in one- 
dimensional coupled cavities. We locate the phase boundary, discuss the behavior of the polariton 
compressibility and visibility fringes across the critical point, and find a non-trivial scaling of the 
phase boundary as a function of the number of atoms inside each cavity. We also predict the 
emergence of a polaritonic glassy phase when the number of atoms fluctuates from cavity to cavity. 



Over the past two decades a considerable understand- 
ing of the physics of strongly interacting systems has been 
gained by a judicious design of controlled many-body sys- 
tems. Successful examples of this sort were optical lat- 
tices or Josephson junction arrays (see the reviews 
The recent proposals 0, 0, 0] to realize a Mott phase of 
polaritons have paved the way to use coupled cavities 
for the study of strongly correlated phenomena in a con- 
trolled way. The rich scenario which emerges in these 
systems stems from the interplay of two effects. Light- 
matter interaction inside the cavity leads to a strong ef- 
fective Kerr nonlinearity between photons. By control- 
ling the atomic level spacings and the photonic resonance 
frequency inside the cavity.it is possible to achieve a pho- 
ton blockade regime [1, 0, @, Q , thus suppressing photon 
fluctuations in each cavity. On the other hand, photon 
hopping between neighboring cavities favours derealiza- 
tion thus competing with photon blockade. Coupled cav- 
ities can be realized in a wide range of physical systems, 
from nanocavities in photonic crystals jlQ] to Cooper pair 
boxes in superconducting resonators [llj . It is therefore 
possible to study a whole new class of strongly interact- 
ing systems that, for the first time, can be addressed and 
measured locally. 

The polariton Mott insulator has been predicted in two 
cases. Hartmann et al. [3J discussed a cavity doped with 
N four-level systems in the limit of large N, while An- 
gelakis et al. [5j and Greentree et al. [4[ analyzed the 
Jaynes-Cummings [l2[ model as a scheme for the light- 
matter interaction; in this last case, an experimental 
proposal has been also devised by Neil Na et al. [HI]. 
Hartmann et al. found a mapping onto a Bosc-Hubbard 
model [IH for the polaritons in the limit of large number 
of atoms and large detuning. In the other case, the phase 
boundary was evaluated at a mean field level for one [J] 
and many [1 31 ] atoms in cavity. The exact phase diagram 
has not been worked out so far; this is what we accom- 
plish in this work for the one-dimensional case. By means 
of numerical simulations and analytical calculations we 
are able to locate the phase boundary and its non-trivial 
scaling as a function of the number of atoms in the cav- 
ity. Furthermore we consider the case where the number 
of atoms fluctuates in each cavity and we show that this 
leads to the existence of a polariton glass. 



The Hamiltonian for the system composed by an ar- 
ray of L identical coupled cavities is given by the local 
Hamiltonian on each cavity and the photon hopping term 
between different cavities: 
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As in 0, [Hj], we add a chemical potential fi. In the previ- 
ous expression t is the nearest-neighbor inter-cavity pho- 
ton hopping, and a, is the photon annihilation operator 
in the i-th cavity; the local contribution 1i!f^ describes 
the light-matter interaction. We will consider the follow- 
ing two models. 

Model I A collection of N two-level systems which 
interact with photons via a Jaynes-Cummings coupling 
H { ^ a) = e(5f + f ) + wajoj + /3(S+ ai + Sr a \), where 
we have defined the spin operators Sf = S^Li G< j,i 
(a = ±, z) and cr^ are the atomic raising/lowering oper- 
ators for the j-th atom, e denotes the transition energy 
between the two atomic levels, u is the resonance fre- 
quency of the cavity, and /3 is the atom-field coupling 
constant (e, w, (3 > 0). The total number of atomic plus 
photonic excitations and the total atomic spin Sf on each 
site are conserved quantities. The ground state is always 
in the subspace of maximum spin, S = N/2. 

Model II In the Jaynes-Cummings model at a large 
detuning A, when the atomic spontaneous emission is 
minimized, also the strength of nonlinearities is weak- 
ened. In order to overcome this problem, a differ- 
ent scheme involving four-level atoms, has been pro- 
posed 15] producing a large Kerr nonlinearity with vir- 
tually no noise. In the interaction picture, in electric 
dipole and rotating wave approximations the model reads 
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SSf+ASt 4 +n(Sf 3 +Sf 2 )+ gi (Sl 3 a^+Sf 1 a) + 
(72 (S? 4 at + Sf 2 a), having defined the global atomic rais- 
ing and lowering operators S Lm — YljLi 10 j ( m \> °j m 
/ (m\ are the atomic raising and lowering operators 
(I to), or energy level populations (I = to) for the 
j-th atom. The transition |3)^ — > |2)^ is driven by a 
classical coupling field with Rabi frequency £1; the cav- 
ity mode of frequency u ca v couples the — > |3) ■ and 
\2) j — > |4}j transitions with coupling constants gi and 
(72; the parameters 8 and A account for the detunings of 
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levels 3 and 4 respectively. The atomic part of the system 
wavefunction for the i-th cavity can be fully character- 
ized by the number of atoms in each of the four possible 
states: n.2, r&3, 714)}, with Ylt=i n i ~ ^ ne total 

number of photons plus the number of atomic excitations 
in the whole system (where states \2) - , |3) . count for one 
excitation, while |4) ■ counts for two excitations), is a con- 
served quantity. Hereafter we assume g\ ~ 52 = 9 and 
define the relative atomic detuning S w = A — S. 

Mott insulator The phase diagram of the coupled cav- 
ity system is characterized by two distinct phases 0,0, @] : 
the Mott Insulator (MI) is surrounded by the Superfluid 
(SF) phase. In the MI polaritons are localized on each 
site, with a uniform density p = n' pol /L, where n po1 is 
the total number of polaritons in a system of L cavities; 
there is a gap in the spectrum, and the compressibility 
k = dp/ dp, vanishes. A finite hopping renormalizes this 
gap, which eventually vanishes at t* . The phase bound- 
aries between the two phases can thus be determined by 
evaluating, as a function of the hopping, the critical val- 
ues of p at which the gap vanishes. Our data have been 
obtained by means of the Density Matrix Renormaliza- 
tion Group (DMRG) algorithm with open boundary con- 
ditions 16]. In numerical calculations, the Hilbert space 
for the on-site Hamiltonian is fixed by a maximum num- 
ber of admitted photons n v J}°* We chose n^b°l = 6 for 
model I and n^°^ = 4 for model 77; we also retained 
up to m = 120 states in the DMRG procedure, such to 
guarantee accurate results, and checked that our data 
are not affected by increasing n^ a °*. We simulated sys- 
tems with up to L = 128, and up to N — 5 atoms per 
cavity 17|; the asymptotic values in the thermodynamic 
limit have been extracted by performing a linear fit in 
1/L. By combining these results with strong coupling 



perturbation theory [18| we were able to locate the phase 
boundaries for all values of N. Most of this paper is de- 
voted to the case e = u> for model I and 6 — A = for 
model II. These regimes could not be accessed by the 
perturbative approach of Ref. [|[ . 

Let us start with zero photon hopping (t = 0). For 
model I, at fixed TV, there exists a value 5} of the de- 
tuning 5 1 = lo — e such that, for Si > <5J, the width 
of the lobe with a polaritonic density p = N is greatly 
enhanced with respect to the other lobes. We estimate 
<5| numerically and find a scaling 6} ~ V7V. For model 
II, at a given relative atomic detuning (5 W > the situa- 
tion is similar to model /, where the resonating lobe with 
p = N is much larger than the other lobes, if 5 < 5*. In 
the opposite case, 5^ < 0, some of the lobes disappear. 

For model /, numerical data at finite photon hopping 
for different values of N are shown in Fig. [TJ the phase 
diagram of model II is shown in Fig. [2] Several interest- 
ing features emerge in the structure of the lobes. In both 
models, for fixed TV, contrary to Bose-Hubbard model, 
the critical values t* of the hopping strength at which 
the various lobes shrink in a point are not proportional 
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FIG. 1: (Color online) Upper panels: Phase diagram for the 
Hamiltonian model I, with AT = 1, 2 atoms inside each cavity 
at e — ui. Lower panels: System compressibility k for the first 
lobe (i.e., p = 1), for different system sizes L, with N = 1 
(left) and N = 2 (right). 
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FIG. 2: (Color online) Upper panels: Phase diagram of model 
//. The detuning parameters are set to zero and g/£l = 1. 
Lower panels: System compressibility k, for the first lobe, for 
different system sizes L, with N = 1 (left) and N = 2 (right). 



to the lobe width at t = 0. Furthermore, the ratio be- 
tween the upper and the lower slopes of the lobes at small 
hopping is greater than the one predicted in Ref. [l8j ]: 
this discrepancy disappears on increasing the number of 
atoms inside the cavity. In terms of an effective Bose- 
Hubbard model, this may be understood as a correlated 
hopping of the polaritons, i.e., the hopping depends on 
the occupation of the cavity. 

A more detailed analysis of the transition at filling one 
can be performed by considering the compressibility. In 



3 



our DMRG simulations we fix the total number n po1 of 
excitations in the system, thus fixing the polariton den- 
sity p inside each cavity in the insulating regime. In the 
lower panels of Figs. [TJ [2] we plot n in the first insulating 
lobe (with a polariton density p = I) as a function oft for 
different sizes of the system, for the models I and II. By 
exploring the mapping to the Bose-Hubbard model we 
construct the full curve t* versus N. The effective repul- 
sive interaction U e s between polaritons, at p = 1, is given 
by U e g(l) = d 2 E(n) / dn 2 \ n= \ (where E(n) is the ground 
state energy of Hamiltonian Hi with n polaritons), that 
is exactly the opening of the first lobe at t = 0. For 
model / it is possible to give an exact analytic formula: 
U eS (l) = 2\/iV[l - y/1 - l/(2N)}, while for model 77 it 
can be evaluated numerically. As N increases, U e g(l) 
decays to zero; for both models U c s(l) ~ l/\/N as far 
as N — ► oo. Moreover the effective repulsion depends on 
the number of polaritons, contrary to the Bose-Hubbard 
model; this dependence weakens, and eventually vanishes 
in the limit n -C N. Therefore the mapping becomes ac- 
curate when N increases. The polaritonic hopping t e g 
can be obtained by performing a strong coupling expan- 



sion in t. For model / we found that t* 
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for model II we get t* 
is then obtained using the value for the critical point 
tlg/UeS — 0.3 jl9|. Figure [3] (upper part) displays both 
numerical (blue squares) and analytical estimates (red 
circles) for the two models. This analysis shows that the 
Bose-Hubbard model provides a good description already 
for N ~ 10. A study of the dynamics is needed to further 
strengthen this observation. We point out that, in exper- 
imental realizations, the parameter that can be changed 
to cross the transition is the detuning. For model I this 
is shown in the lower part of Fig. [3l 

Visibility of photon interference The phase transi- 
tion can be detected by analyzing the phase coherence 
of photons [2l| . in a way similar to what has been 
done 



21|, I23. The 



for the Bose-Hubbard model 
terference pattern of the photonic density is propor- 
tional to the photon number distribution S in the mo- 
mentum space: S(k) = ± T,j,i=i e 2mk ^' l ^ L (oto,) . The 
visibility of interference fringes can then be defined as 

V = (S max -S min )/ (Smax+Smin) • The visibility is strictly 
zero only in the limit t = 0, where the interference pat- 
tern S is constant. When t is increased, the visibility 
itself increases, until it saturates to the maximum value 

V = 1 in the superfluid regime. This description is well 
adapted even for photon coherence in our hybrid light- 
matter system, as shown in Fig.|U The existence of differ- 
ent phases can also be detected by measuring fluctuations 
in the number of polaritons, as discussed in 

Polaritonic glass phase Up to now we have assumed 
that the number of atoms in each cavity was constant 
and equal to N. In certain implementations this require- 
ment might be demanding. Here however we consider 
this problem from a different perspective and show that, 
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FIG. 3: (Color online) Upper graphs: Critical hopping t* in 
the first lobe as a function of the number N of atoms inside 
each cavity, for the two different models with all the detunings 
set to zero. Dashed black lines indicate a behavior t* ~ 1/ \fN 
and are plotted as guidelines. Lower graph: Critical hopping 
t* for model / with N = 10 atoms per cavity [201 ] . as a function 
of the relative detuning A. Blue squares have been evaluated 
with the DMRG. Red data are estimates obtained from the 
effective on-site polaritonic repulsion U e a at zero hopping, 
with filling p — 1. 



when N changes from cavity to cavity, it leads to the 
emergence of a polariton glass. Following , this phase 
is characterized by a finite compressibility, gapless exci- 
tation spectrum, and zero superfluid density. 

Random fluctuations in the number of atoms per cav- 
ity lead to disorder in the on-site light-matter interaction 
strength. This effect can cause significant consequences 
only in the limit of large N (we quantify this state- 
ment below) where the mapping onto a Bose-Hubbard 
model applies. A Bose glass phase has been originally 
predicted [l4[ as a function of disorder in the chemical 
potential. Recently it has been shown that fluctuations 
in the on-site repulsion lead to a Bose glass as well j23|. 

We take advantage of the results obtained in and 
the mapping to an effective Bose-Hubbard model to give 
a detailed estimate for the width of the polariton glass 
phase. The key to our finding is to relate fluctuations 
in the number of atoms to fluctuations in the on-site 
repulsion U c ff. We suppose that each AT$ is a random 
discrete Gaussian variable with a mean value (N), and 
a standard deviation SN. Figure displays fluctuations 
in the effective on-site repulsion 5U e g as a function of 
the fluctuations in the number of atoms, while in the 
inset of Fig. [5] we show an example of such a varia- 
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FIG. 4: (Color online) Photonic visibility V in the first insu- 
lating lobe (p — 1) as a function of the hopping strength tj '13, 
for model I. Here we considered the case JV = 1, A = 0. In- 
sets: The photon number correlation function S(k) for L — 48 
in the deep insulating regime (t/(3 = 0.01), and at the edge 
of the superfluid regime (t//3 = 0.20). 
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FIG. 5: (Color online) Relative standard deviation of the ef- 
fective on-site interaction strength U e g averaged over L = 10 4 
cavities, as a function of the atom number fluctuations SN; 
from left to right (N) = 1, 10, 100, 1000. Dashed line indi- 
cates an effective variation of the interaction strength equal to 
the standard deviation of a random interaction uniformly dis- 
tributed in the interval U c s(i)/ (U e g) 6 [— e, +e] with e = 0.25. 
In the inset an example of the variation on the on-site effective 
interaction is shown; (N) = 100, SN = 20. 



tion, with (JV) = 100, SN = 20, in the case of model 
/. We can then use numerical data of Ref. [23| to esti- 
mate where the polaritonic glass phase can be observed. 
It has been shown that in the Bose Hubbard model, for 
a relative interaction 5U c s/ (U c n) uniformly distributed 
in an interval of length 2e = 0.5, a L = 200 sites sys- 
tem at filling p = 1.01 exhibits a Bose glass phase for 
0.078 < t e f{/ (U e f{) < 0.133. In our model with the same 
polaritonic filling, if we take, e.g., (JV) = 100 particles per 



cavity and we choose SN ~ 19 (such to have e = 0.25), 
a polaritonic glassy phase should be visible, in the case 
of model I, for 7.51 x 10~ 3 < t/f3 < 1.401 x 10~ 2 . The 
situation is qualitatively the same for model II. The 
polariton glass may be observed by measuring, at frac- 
tional fillings, the photon visibility as a function of the 
disorder. Anyway, a complete characterization of the po- 
lariton glass requires a study of its dynamical behavior. 

Concluding remarks In this work we have discussed 
in details the equilibrium properties of a chain of cou- 
pled cavities. The results we presented here are general, 
in the sense that they apply to all the systems described 
by Eq. |T]). An important issue is to understand the ef- 
fect of decoherence, decay, cavity losses, . . . , which occur 
in these systems. These aspects have to be discussed for 
each implementation (see for example 0, [ll[ ) , together 
with the effect of an additional possible external driv- 
ing. The combined presence of dissipation and external 
driving may lead to a change of the universality class of 
the transition, or to new phenomena associated to non- 
equilibrium phases. 

We acknowledge useful discussions with A. Tredicucci 
and support by MIUR-PRIN and EC-Eurosqip. 
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